What does the oscillate into? 



Debajyoti Choudhury 1,2 and Anindya Datta 2 
1 Department of Physics and Astrophysics, University of Delhi, Delhi 110007, India 
2 Harish- Chandra Research Institute, Chhatnag Road, Jhusi, Allahabad 211019, India. 

The favoured resolution of the atmospheric neutrino anomaly involves an oscillation of the muon 
neutrino to a different state. Current experiments allow for the latter to contain a significantly large 
fraction of a non-standard flavour. We demonstrate how the next generation of experiments may 
take advantage of matter effects to resolve this issue. 
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A combination of several seminal experiments [1|, 

have left us with the inescapable conclusion that 
neutrinos oscillate amongst themselves. We do know 
that the muon-neutrino (z/^) mixes almost maximally 
(sin 2 26» m3 > 0.92 at 90% C.L) with another species vz 
(distinct from the electron- neutrino u e ) and that the 
mass splitting |<5m 2 2 | ~ 2 x 10~ 3 eV 2 . Furthermore, u e 
mixes with a combination of and v T with a similarly 
large angle 9 e 2 ~ 30° but a far smaller mass splitting 
(5m 2 ! ~7x 10~ 5 eV 2 ). The remaining mixing angle #13 
is constrained to be very small |3j . It is often touted that 
neutrino oscillation physics is now poised to move into 
the domain of precision studies, with the major remain- 
ing questions primarily centering on unraveling the sign 
of <5m| 2 , accurate measurement of all the parameters and 
looking for possible CP violation in this sector. 

Several other questions remain though, not the least 
of which relates to the identity of v 3 . For one, this could 
be a linear combination of v r and a sterile neutrino (ys)'- 



f 3 = C£ 



S£ Us 



where = cos£ and sj = sin£. As long as 6 3 = 6^3 ~ 
7r/4, the above choice would reproduce the charged cur- 
rent (CC) data at both Super-K Q and SNO % While 
data already rules out S£ = 1, note that s 2 = 0.25 (0.35) 
is still allowed at 90%(99%) C.L. @. In fact, v s could 
even be an antineutrino (say, the v T ) Q if one would 
allow for a violation of the total lepton number. If, in 
addition, a helicity flip is enhanced (being proportional 
to m v , this is normally suppressed) or a new interaction is 
relevant, this could even produce an active v r , and thus 
induce wrong-sign CC events (analogous oscillations of 
the form v e ,iAj, — > v e have been looked for at the BEBC 
experiment |8| with no evidence being found). Unfortu- 
nately, none of the current or forthcoming experiments, 
including OPERA [j| and ICARUS E3i are well-equipped 
for such measurements. 

In this article, we seek to point out how passage 
through dense matter may amplify the effects of a non- 
zero so that CC events can be used as a discriminator. 
Since the refractive indices for different neutrino beams 
may be altered to a significant degree E3 by their in- 
teractions with matter, clearly the survival probability 



Pup, will display non-trivial matter effects for <-> v$ 
oscillations while it would not do so for <-> v T oscil- 



lations. For an arbitrary s^, we will see that the resul- 
tant P[i[i would not necessarily lie in between the two 
extreme cases. Rather, the latter has a non-trivial de- 
pendence on sj. To quantify this, let us begin by re- 
stricting ourselves to the relevant three-neutrino sub- 



space namely iv j 



two mass (square)-differences Ai^. 
tonian can then be expressed as 



vs), characterized, in general, by 



H 




The effective Hamil- 
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where Ji is symmetric and Si = Ai/(2E V ). Vno — 
—Gf N n /2^/2, with N n denoting the neutron density in 
matter, denotes the relevant matter effect. 

Note that, for a non-zero sj, matter effect arises even 
in the absence of any mixing with the v e . On the other 
hand, for = 0, all matter-induced effects (including 
(1) those in the v tli —v T sector) are proportional to the size 



of the relevant mixing with v e and are understandably 
small 

[3 El 111 13 El 13 El The approximation of 
a three-generation mixing scenario is, thus, an excellent 
one and has the virtue of exhibiting the essence of the 
effect. 

Although the effective mixing angle now becomes a 
function of instantaneous density, it should be realized 
that it is not necessarily enhanced on account of the mat- 
ter effect. Rather, for sj = 1, when the v T decouples, the 
Vfi—fS mixing, in vacuo, is close to maximal and the mat- 
ter effect only serves to decrease its size. On the other 
hand, for small S£, a resonance enhancement of the v^-vs 
mixing is indeed possible. In Fig. ^ we exhibit the sur- 
vival probability for atmospheric of a fixed energy as a 
function of its incident angle. For the numerical results, 
we have chosen #23 = 7r/4 and used the varying density 
profile of the earth as give n by the Preliminary Reference 
Earth Model (PREM) [ijj . The aforementioned suppres- 
sion, for Sf = 1, is clearly evinced by the corresponding 
P w never vanishing identically even at the minima (in 
stark contrast to = 0). For intermediate S£, the two- 
generation simplification does not work, and the change 
in P^ is more subtle. For small E v , where the diameter 
of the earth permits many oscillation periods, the value 
of Pita at the maxima decrease monotonically. While this 
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-0.9 -O.V -0.5 -0.3 -O.l FIG. 3: (a) The number of f M CC events, as a function of 

cos e incidence angle, from atmospheric neutrinos (with 5 GeV < 

FIG. 1: The atmospheric survival probability as a function E v < 40 GeV) in an ICAL m/ like detector. The mixing 

of the incident angle. Each panel refers to a different E v . The angle # M 3 = 7r/4 and Ai = A2 = —0.003 eV 2 . (b) The 

mixing angle 6^3 — tv/4 and Ai = A2 = —0.003 eV 2 . weighted asymmetry between and CC rates. 



effect scafes with for small £, for large sg, the system 
quickly becomes rather nonlinear. For large E v , the mat- 
ter effect is understandably much larger and non-linearity 
sets off even for small . 

A detector-independent measure of the sensitivity 
is afforded by the quantity 0] 1Z{E V , cos 0, s%) = 

KM - Pw(°)] I [ p »M + p w(0)] 1/2 - The contours 
of constant TZ (FigEJ, apart of illustrating the com- 
plicated dependence on also indicate the choices of 
(E v , cos 9) ideal to probe a given . Note that high val- 
ues of 1 72. 1 are also obtainable at E v < 5 GeV, but owing 
to the rapid oscillations in P w , the exclusive use of such 
neutrinos would require very fine detector resolutions. 

Such effects are characteristic of the scenario and can 
thus be used even to measure S£. Furthermore, it should 




FIG. 2: Contours of constant TZ (defined in the text) in the 
E,j — S(. plane. Parameters are as in FiaYll 



be realised that with Vnc being of the opposite sign for 
antincutrinos, this deviation from the simple = case 
would occur for only one of neutrinos and antincutri- 
nos. Thus, if the detector were a magnetized one — such 
as the proposed iron calorimeter (ICAL) at the India- 
based Neutrino Observatory (INO) [2(J — and hence able 
to measure the muon-charge, this could be used as a fur- 
ther check if any such deviation were to be seen. Of 
course, whether z/^'s show this effect or u^s would de- 
pend on the sign of the mass difference and thus would 
constitute a measure of the latter were to be non- 
zero. However, note that at large E v , where the differ- 
ence is more pronounced, the atmospheric neutrino flux 
falls rapidly. On the other hand, at lower energies, one 
would need both very good energy resolution as well as 
angular resolution to draw any positive inference. 

As an example, we consider the abovementioned ICAL, 
where the projected energy threshold for muon detection 
is about 2 /gev and the resolution is expected to be better 
than 0.5 /gev over the entire range |18L l20l | - In Fig|3i we 
show the number of \x~ events that such a detector would 
accrue from atmospheric neutrinos |2lj | for a 300 kT-yr 
running period. To account for the threshold, as also the 
fact that the difference is more apparent for larger ener- 
gies, we integrate over the range E v > 5 GeV. Clearly, 
the present limit of < 0.5 would be easily distinguish- 
able from the — case. In fact, a sterile component 
one-fourth as small (s^ = 0.25) would begin to be visible. 
Note that the v events are not expected to show a similar 
effect for negative A,. It is thus interesting to consider 
an asymmetry between the ^ event rates N± 

A K = (N- -kN+)/(N- + k N+) . 

Using k ~ 2.09 to compensate for the difference in 
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and CC cross-sections over the entire range, the cor- 
responding A K is displayed in Fig©. Clearly, this too is 
a good discriminant, though not as powerful as AL. itself, 
owing to the smallness of N+. 

Although the resolving power of an ICAL-like detec- 
tor using atmospheric neutrinos has been demonstrated 
above, it is still contingent upon us to identify possible 
experiments with a greater sensitivity, in particular to 
probe the case of positive Aj. Thus, it would be par- 
ticularly useful if one could avail of a beam of high en- 
ergy neutrinos. The latter is possible in the context of 
accelerator experiments, especially a future facility such 
as a super-beam or a neutrino factory wherein an in- 
tense beam of muons is to be accelerated to a not too 
high energy and stored in a storage ring with a straight 
section directed towards a neutrino detector 22]. With 
muons decaying in this straight section thereby produc- 
ing high intensity neutrino beams (both of electron- and 
muon-types) that are highly collimated in the direction 
of the decaying muons, such facilities have two additional 
advantages over neutrino experiments with conventional 
neutrino beams (arising from ir^ decay), namely (i) a 
precise knowledge of the and v e fluxes helps reduce 
the systematic errors; (ii) assuming a [i-beam of, say, 
20 /gev energy, such neutrinos have, on the average, en- 
ergies higher than those of conventional neutrino beams 
thereby increasing the cross-sections at the detector. 




FIG. 4: The survival probability as a function of E u for 
different baselines. Parameters are as in i*%Ql 

In Fig. 01 we display the survival probability for such 
accelerator neutrinos as a function of their energy for a 
given baseline. As expected, the difference between the 
two scenarios grows with the baseline, for it is only with 
large baselines that the neutrinos sample the larger den- 
sities in the earth's interior. With the differences starting 
to become apparent only for baselines larger than 2000 
km, a facility such as ICAL/INO [2(j would be ideally 
placed to make such measurements particularly with neu- 



trinos from JHF (baseline - 5000 Km), CERN (~ 7000 
km), or Fermilab (~ 10000 Km). For a given baseline, the 
difference between the scenarios is expectedly more pro- 
nounced at larger neutrino energies. Furthermore, larger 
E v translates to larger cross-sections and, hence, larger 
event numbers. On the other hand, the dependence of 
the neutrino flux on E v needs to be taken into account. 
Convoluting all such effects, we show, in Fig0 the ex- 
pected number of events for a typical iron calorimeter 
detector [2(j. Although a larger is preferable, both 
on account of larger cross-sections as well as a tighter 
beam collimation, we choose to work with a likely first- 
generation configuration, namely E^ = 20 GeV and a 
total of 10 21 decaying muons. 
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FIG. 5: The number of events expected for a 20 GeV muon 
storage ring and a 50 kT iron calorimeter with a energy 
threshold of 2 GeV for the ^ '2C]. 

Working with the assumption of e s = 0, we neglect, in 
FigEl any wrong sign muon event (which, anyway, can 
be distinguished by a magnetized detector) caused by 
the v e from pT decay oscillating into while travers- 
ing the earth. For a baseline of 5000 km, the sensitiv- 
ity is marginal and = 0.5 just about separable from 
S£ = 0. At 7000 km, on the other hand, one can eas- 
ily explore down to sin£ = 0.25. For 10000 km, things 
improve dramatically and even sin£ = 0.25 is remark- 
ably distinguishable with just one-tenth the number of 
muons (FigEJ. A further countercheck is afforded by the 
fact that, for Si < 0, one would expect a matter effect 
only for neutrinos (i.e., \i~ run) and not antineutrinos 
(a /i + run). However, due to the lower cross sections for 
v^, achieving similar significance in the negative result 
would require twice the running time. 

It should be realized that some of the details of 
Figs (TjyU are a consequence of our choice of 0^ = ir/4. 
As FigElshows, were this mixing to be non-maximal, 
would be strictly positive for the canonical case of = 0. 
On the other hand, for Sf = 1, at resonance, the effec- 
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FIG. 6: The neutrino survival probability as a function of the 
neutrino energy for different baselines and sin 2 2# M 3 = 0.92. 
The parentheses in the legend refer to (s^, Aj) combinations. 
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FIG. 7: As in Fig\j% ftitf /or sin 2 26»^ 3 = 0.92 instead. 



tive v^-vs mixing would get enhanced, thereby allowing 
Pfifi to reach zero at the minima. For intermediate sj, 
the two effects get intermixed leading to a damping in 



the oscillation at both ends (minima as well as max- 
ima). Note that a non-maximal value of 6^3 renders the 
sign of Si very important. Were 5i to be positive instead, 
clearly there would be no resonance enhancement of the 
transition probability even for S£ = 1. 

The corresponding event expectations are displayed in 
FigEI Clearly, deviation of 6^3 from maximality ren- 
ders a non-zero value of even more visible, at least for 
Aj < 0. A positive value for the latter understandably 
reduces the significance to the extent that, even with 
10 21 muons, the case of = 0.25 is just barely dis- 
cernible. However, it should be kept in mind that, in 
such a case, antineutrino oscillations have an enhanced 
sensitivity and the effect would be clearly visible even 
accounting for a reduced cross section. 

In summary, we have exhibited that matter effects may 
be used profitably in resolving the vexing issue of the 
identity of the neutrino state that oscillates into, es- 
pecially since the former is still allowed to contain a sig- 
nificantly large fraction of a sterile state (or, equivalently, 
an antineutrino state). Studies with atmospheric neutri- 
nos may themselves be used to obtain a reasonably high 
degree of resolution. For example, with just 300 kT- 
yr exposure, the planned ICAL/INO detector should be 
sensitive to a sterile admixture one-fourth the size of cur- 
rent bounds. With the use of a muon storage ring with 
a very long baseline, the resolving power improves dra- 
matically. Furthermore, such experiments are shown to 
be quite sensitive to the deviation of the effective atmo- 
spheric neutrino mixing angle from maximality thereby 
allowing for an accurate measurement of the same. And 
while we have neglected the mixing with v e in our anal- 
ysis, such effects have been demonstrated to be small 
and would not lead to any significant deterioration in 
the analysing power. 
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